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Abstract 

The extended principle of minimal action is described in the presence of prescribed 
, source and sink points. Under the assumption of zero net flux, it leads to an optimal 

Monge-Kantorovich transport problem of metric type. We concentrate on action corre- 
sponding to a mecahnical Lagrangian. The optimal solution turns out to be a measure 
p ^ ' supprted on a graph composed of geodesic arcs connecting pairs of sources and sinks. 

^ ' 1 Introduction 

1.1 The extended action principle 

Consider the Lagrangian of a mechanical system: 

Lip,x) = hp\^-V{x) , 

■ where x,p G M and V is a smooth bounded potential function. 

. The minimal action is a function on M'= X M'^' X M X M, defined by 

in 

J{xo,xi,to,ti) :=mi I L{y{t),y{t))dt , (1.1) 

where the infimum is taken on all orbits y : [to,ti] — > satisfying the end conditions 

yito) = xo, y{ti) = XX. 

The extended minimal action principle (EMAP) was introduced by Benamou and Brenier 
([BB,BBG]), in the case 1/ = 0, and by [W] in the general case. Below we review its definition 
' and relation to the differential point of view of optimal mass transportation - see also ch. 8 

of [V]. _ 

Let M be the set of probability Borel measures on R*^. Let £ the set of M*^— valued 
Borel measures defined on R'^. The extended Lagrangian is defined for any {^x, E) ^ M. x £ , 
provided E is absolutely continuous with respect to /i and the Radon-Nykodym derivative 
dE /dn G L^((i/i), as 

'dE 



X 



^E,li) ■■=1 



dfi 



Vix) G Li(d;u). 



The extended action is defined on all orbits fj, = fit (res. E = Et) oi J\A (res. E) valued 
functions of the real line, which satisfies the continuity equation 

dtfit + • = (1.2) 
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in the sense of distributions. For any pair //Q) /^i & -M. and any to; G M we define tlie EMAP 
J{fio, fxi,to,ti) := min / / C ( E t, fj-t) fJ-t{dx)dt (1.3) 

{^l,E)Jto JM*= ^ ^ 

over all pairs ^/Ut, Ei^ satisfying (|1.2() . subjected to the end condition /if^ = /io, /Uj^ = /Ui. 

EMAP is really an extension of the classical minimal action principle, if the end measures 
^Oi/^i are replaced by the point measures S(^x^y In general, a minimizer {Et,^J-t) of 

the EMAP exists and the Radon-Nykodym derivative dEt/dfit is a Lipschitz function on 
Supp{fi) under some general conditions, so that the flow induced by 



dEt 



dfit 



(1.4) 

xit) 



exists and is unique for fi a.e. point {x,t). Moreover, equation H1.4() is compatible with the 
Euler-Lagrange equation associated with the Lagrangian L for each individual orbit x = x{t) 
(see [W]). 

The motivation of Brenier and Benamou for the introduction of the extended action prin- 
ciple is an attempt to devise an algorithm for solving the Monge optimal mass transportation 
for quadratic cost: Given two probability measures fio^fJ-i on a common space (say R'^), find 
a mapping T : R*^ ^ R*^ which transports fio to /ii so that the cost of transportation 

min / |T(x) — xpd/iQ • (1-5) 

Recall that a mapping T transports //q to fii {T^^q = /ii) if and only if (^~^(^)) = ^i(^) 
for any (Borel) measurable set A. The existence of a unique minimizer of (jl.Sj) is known 
if /io is absolutely continuous with respect to Lebesgue measure and both /iOjA*! have finite 
second moment (see [B] , and later extension in [GM] ) . It was the fundamental observation of 
Benamou and Brenier that (in the case V = 0) the flow (|1.4|) associated with the solution of 
the EMAP induces the optimal mapping T of H1.5|) under some regularity assumptions. The 
extension of this result to the case y ^ 0, introduced in [W], relates the flow (|1.4|) to the 
optimal solution of the mass transport with respect to the cost function c = c{x,y), 

min / c{x,T{x))dfiQ (1-6) 

where c{x,y) = J{x,y,ti,t2) given by the classical action 

1.2 Sources and sinks 

In this paper we consider the extended action principle under a set of sources and sinks. Let 
Di C R'^, i = l,...n be a disjoint set of compact domains whose boundaries Si = dDi are 
smooth surfaces. On each such surface we assign an integrable function Aj : — > R so that 



N 

L 

=1 



(p Xids = . (1.7) 
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The extended action principle under the presence of prescribed fluxes Aj across the surfaces 
Si is defined as H1.3|l . where the set of pairs ^/i, ^) is defined on (M'^ - U^A) x [to,ti] and 
(|1.2|) replaced by 

dtfi + V^-E=0 on (m'^-U^a) x[to,ti] , Efni = Xi; V(s, t) e 5^ x [to, ti] , (1.8) 



where nj(s) are the outward normal to Si at point s. 

We shall concentrate on stationary extended minimal action: nt ( and, correspondingly, 
Et), are independent of t. 

Under stationarity condition, the extended principle of minimal action takes the form 



mm 



-UD,; 



dE 



dfj. 



Vdfi 



(1.9) 



where £ M is supported on M'^ — UDj and E is a vector valued measure, absolutely 
continuous with respect to /i, which satisfies in addition 



V • ^ = on 



E ■ rii = Xi on Si, 1 < i < n . 



(1.10) 



Points sources and sinks: Letting the domains Dj shrink to points xj € M*^, we replace 
the functions Aj on Si by constants Aj € M so that 



The admissibility condition ()1.1U|1 for the stationary minimal action H1.9() can now be casted 
into the single condition: 



V-E=^Xi5^^ 



on 



(1.11) 



The problem is now reduced to finding the minimizer G of 

2 



:= min 



dE 



dfj. 



dfi 



Vdfi 



where E is absolutely continuous with respect to and satisfies p.lljl . 



:i.i2) 



1.3 Main results 

In this paper we concentrate on point sources and sinks. The set of fluxes Ai, . . . A„ is divided 
into two sets: sources (Aj > 0) and sinks (Aj < 0). 
Let 

/+ := {i;Xi>0} ; /_ := {i ; A^ < 0} . (1.13) 
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The metric Monge transport plant between a pair of probability measures HojfJ-i is the 
special case of (|1.6() where c{x,y) is a metric, say c{x,y) = \x — y\. It turns out that the 
solution of the stationary EMAP is related to the metric Monge transport where /xo,^i are 
replaced by the point measures 

where |A| := = - ^i- 

In general, however, there are no mappings T which transport a point measure /xq to 
another measure, so there is no sense to define the Monge problem 1)1.6^ for this case. However, 
the Monge problem can be relaxed to an optimization problem on the set of 2-point probability 
distributions Q = Q{dxdy). This is the celebrated Kantorovich relaxation of the Monge 
problem [K]. 

We shall first describe the special case V = 0, corresponding to the Euclidian metric 
c{x,y) = \x — y\. In the case of point measures (|1.14() . the Kantorovich relaxation of the 
minimal cost takes the form of the 1— Wasserstein metric (see, e.g. [R]). It is defined by 
a |/+| X |/_| matrix A° which realizes the minimum 

w^'H^) = E E - ■= E E ^mI^^ - ^ii • (1-15) 

where 

•= I > ; E = Ai i G /+ ; E ^^.i = \ ■ (1-16) 

The first result states that the minimal action 1)1. 12() is given by 
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Let Go be the bi-graph composed of the set of vertices Xi, i = 1, . . . n, and edges L{i,j) 
are segments connecting i G /+ to j G /_. 

The second result states that the action minimizer of 11.12)1 is supported in the 

edges of the bi-graph Go for which A" j > 0. On each such edge, ^ is a uniform measure and 

^^{Li^,j)) = \X\-'Al^. 
We now consider the case V ^ 0. We assume that 

lim V{x) = . (1.17) 

\x\—too 

Let 

V = supV . (1.18) 



4 



For any E > V we consider the Riemannian metric associated with the Maupertuis' action 
principle (see [Ar]): 

daE = VE - Vds , (1.19) 

here ds is the Euclidean metric. 

The geodesic distance associated with this metric is denoted by De{-, •)• We recall that 
a geodesic arc connecting two point x, y coincides with an orbit of the mechanical system 

X + VxV{x{t)) = 2;(0) = xo, x{T) = xi , 

corresponding to the energy level |xp/2 + V{x) = E. Here T = Tx^^xi (E) is the time of flight 
from Xo to xi (which is, of course, a function of E as well). 

The bi-graph is defined, analogously to Go, as the collection of vertices Xj, 1 < i < n, 
and all edges composed of geodesic arcs (with respect to the metric das) connecting Xj, 
2 G /+, to Xj, j £ I . The 1— Wasserstein metric associated with this distance is given, 
analogously to W'^-'^^ as 

Wg^(^) = E E Af^jDEix^,x,) := ^min ^ ^ ^„Z)s(x„ x,) . (1.20) 
Then the minimal action 1)1.12(1 is 

' ^ ' (1.21) 



= max 

E>V 



V2W« ( A ) - E 



There exists a minimizer /x G realizing this action which satisfies the following: 

Case a: Eq > V is the maximizer of the RHS of (|1.21|) . Then there exists an action 



minimizer supported on the bi-graph Ge so that fi {LE{i,j)) > only if Af.j > 0. To wit: 

lEorp 



^(Ls„(i,j))=2-i/2^f°r,,,(i?o) , (1.22) 
where Tij{Eo) = Tx^.x^iEo) as defined below ()1.19|) . 



Case b: li Eq = V is the maximizer of ([1.2111 then the following holds: Let (j,q be the measure 
supported on Gy subjected to (|1.22ll . Then there exists (3 G [0, 1) so that 

Let Xo be a maximizer of V , that is, V = ^(xo). Then 

/i = /io + 13^X0 

is an action minimizer. 

In both cases (a) and (b), the following claim is valid: 

Time/Flux duality: The expectation of the inverse flow time, (T^^), is proportional to 
the total in(out) flux |A| := — ~ ^j- 

E^(T-i) := E E f^{LEo{i,j))TrliE,) = ^ . (1.23) 
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1.4 Outline 



In section 121 we derive the weak formulation of the stationary EMAP, which leads to a dual 
problem: 



inf sup 



'-[ 









^{dx) + w 4'{s)Xi{s)ds 



In section 131 we concentrate in the case of point sources and sinks, where Di shrink to points 
Xi. It contains some definitions and a preliminary lemma on the Wasserstein metric and its 
dual representation. 

Section jll Is the most technical part of this paper. It contains a sequence of auxiliary 
lemmas, which are needed to the proof of the main result, as described in section 1131 The 
proof itself is given at the end of this section. It is given for the general case y ^ 0, since 
the case V = Q follows easily form the general one. Finally, the short section summarizes 
the results of this paper. 



2 Weak formulation 

2.1 Notations 

Given n disjoint compact sets Di C and Si = dDi smooth surfaces, set 



Let M stands for the set of probability Borel measures on J7. Let 
M := Iv ; v{dxdp) Borel Probability measure on 1] x M'^ , / 



\p\'^u{dxdp) < CO . 



Let P : J7 X M'^ — > Q be the projection F{x,p) = x. For v ^ M., P#i^ is the push-forward to 
the marginal measure // = fj,{dx) £ M: 



Let also 



For ij, e M let 



'>l){x)v{dxdp) = / i){x)¥#v{dx) ; M^j G C{VL x 



pv{dpdx) 



A := 



IveM; / 'Ey{dx) ■ V0(x) - V </ 



E^{dx) :-- 

s)Xi{s)dpds = V(/) G C\W') 



Af, := An Mf, . 
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The Lagrangian is now defined as a function on M. via: 



L{u) = ]- j \p\^u{dxdp) - [ V{x)F#{u){dx) 
2 7nxM'= Jn 



2.2 Weak form of the minimal action 

We now describe the weak form of the stationary action principle: 

minL(z/) . 
ueA 



(2.1) 



(P) 



Equivalently, if 

then (P) is equivalent to 



L(/x) := inf L{u) 



'mj_L{iJ,) . 



2.3 Dual representation 

Let now define C : M x C°°{n) ^ RU {oo} by 



£{u,(j)) := L{i^) - f E^{dx) ■ V(t>{x) + V / (j)Xi ■ 
Jn ^ J Si 

Next, we use an appropriate version of the minmax principle to obtain the dual formula- 



tion: 



Proposition 2.1. For any ji G M, 



Proof. First, note that 



L(/x) = sup inf C{iy,<p). 



L{ijl) = inf sup C(y, > 



0eci(n) 

Indeed, by definition, sup^^^;,!^^ Ciy, 0) = oo if i/ ^ A, while Ciy, cp) = L{i') otherwise. Next, 
note that L is an affine function on each of the domains Ai^ and C^{Q,), separately. As such, 
it is a convex functional on A^^ and concave on C^{Q.). In addition, M.^ is a compace set 
with respect to the weak topology (in which C is continuous). The Minmax theorem, then, 
can be applied (see, e.g. [Ro]), and the claim follows. □ 



Next we evaluate 



inf C{i', (f>) = inf 



1, 



inf 



iy{dxdp) + 2. f 4>{s)Xi{s)dpds . 
u{dxdp) + (p <p{s)Xi{s)dpds . 

A J Si 
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(j){s)Xi{s)dpds . 



(2.2) 



where the infimum is obtained at v{dxdp) = f^{dx)5(^p_-^fj,y 
Let us now define, for any fi £ Q and (p £ C^(r2): 



fj.{dx) + 



{s)Xi{s)ds . 



By Proposition 12. II we get 

Corollary 2.1. L{fi) = sup^g^i^j^) J^{4>) 

Proposition 2.2. // ^1 is compact, then there exists /i G 7W which solves problem P*. 

Proof. If i7 is compact, so is the set Ai with respect to the weak topology, as the set of 
Probabihty Borel measures on a compact set. In addition, L is lower semi continuous, since 
it is a supremum of the affine functionals J[.){4>) by Corollary 12.11 Hence, a minimizing 
sequence of L in contains a subsequence which converges to a minimum of L. □ 



3 Point sources and sinks 

Assume now that the surfaces Si degenerate to points Xi S R*^. The fluxes functions Aj 
defined on Si degenerate, then, to constants Aj G M. The total flux condition takes the form 

n 

Y.Xi = 0- (3.1) 

1 

In this case, the functional is written as 



JlRfe 2 ^ 



[Xi 



Recall that be the distance metric induced by the Riemannian metric daE H1.19() . 
that is, for x,x M^: 



De{x,x'):=M JE-V{q{s))ds (3.2) 
<? Jo 



where the infimum above is taken over all orbits q{s) in Euclidian arc-length parameterization 
connecting x io x . A minimizer orbit q of (|3.2|) is called an E— geodesic for the metric 
De{x,x). _ 

Given E > V (|1.18j) . i,j G (1, . . .n), let qij{s) be the i?— geodesic of DE{xi, Xj). parame- 
terized by the Euclidian arc- length. So, qij{0) = Xi, qij{Sij) = Xj where Sij is the Euclidian 
length of the DE{xi,Xj) geodesic curve. Let 

T^AE) = r'\E-V{qi,{s))r^'^ds . (3.3) 



Note that Tij{E) is the time interval of existence of the orbit which connects Xj to xj at 
energy E. By differentiation of DE{xi,Xj) with respect to E we obtain 
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Lemma 3.1. If E > V then 



-^DE{xi,Xj) = ^Tij{E) < oo . 



For each E > V define 



P^As) ■■= -——j===== ■,0<s< Si, . (3.4) 

Then pij is a probability density on the interval [O.S'ij]. Then 

Atjj := iPi,j{s)ds) (3.5) 

is a push forward of this probability density to a probability measure Hij supported on the 
geodesic arc connecting Xi and Xj in M.^. It is defined, for any test function cf) E C(M'^), 

via 

r rSi,3 

/ (j){x)dfiij = / (p{qi,j{s)) pij{s)ds . 
Jr'' Jo 

We now recall the definition of the 1— Wasserstein metric (jl.2U|) . By duality formulation 
of the Wasserstein metric we also obtain (see [R]): 

W«(A):= max |f;A,<A. ; ^^^^^ < 1 , V j, A; e{l,...n}.\ . (3.6) 

^gRn DE[Xj,Xk) J 

Below we collect some useful results: 

Lemma 3.2. A & Q-^ HI. lb]) is an optimal solution of \1.2U\) iff there exists a minimizer cf) 
of itu/. 6|) so that Aij = for any pair i,j satisfying — < DE^Xi, xj). 

Proof. For any (p satisfying maxj^j — DE{xi,Xj) < 1 and A S we obtain 

XI XI ^i,jDE{Xi,Xj) > X X ^i'Mi ~ ^j) 

ie/+ \je/+- / iG/_ \iG/+ / iG/+ is/- i=l 

where we used Aij = if either i,j £ /_(_ or i,j £ Since equality holds for the optimal A 
and (j) , it follows that Aij > implies (pi — (pj = DE{xi,Xj). □ 

Recall the definition of the subgradient of a function /i at x G M*': 
d^h = ^peR^ ■ h{x') - h{x) >p-{x'-x) Vx' eR^ 
The following result can be found in, e.g., [HL]: 
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Lemma 3.3. 

i) If h is a convex function defined on M'^, then dxh for all x G M'^. 

a) If h is convex, then p G dxh if and only if x £ dph* where h* is the Legendre transform 
ofh. 

Hi) If h is convex, then for any x,p G M'^, the inequality h{x) + h*{p) > x ■ p holds with 
equality if and only if p £ dxh if and only if x £ dph*. 



4 Proof of the main result 

We shall prove the main result of section 11.31 for the general case V ^ satisfying H1.17j) . 
The special case V = follows easily. 
Let ~$ := (t>n) G and 



B- 



For fj.£ M set 



and 



H{ (j) ; fj,) := inf 



I [ \V<P\^dfi+ [ Vdfi 



Hic^) := su_p . 

Remark: Note that H{-;fi) is not necessarily a convex function on M.^ for each fi £ A4. 
However, H is convex, as we shall see later on in Corollarv 14.41 

Lemma 4.1. For any G M", 

= mini ^, max '^^^ < ^2} . (4.1) 

E>v [ DE[Xi,Xj) J 

Proof. Let E > V. Assume there exists (p £ C^{M.^) so that 

|V0| < V2VE - V on and </>(xi) = , l<i<n . (4.2) 

Now, if such a function satisfies 1)4. 2|1 . then for any i ^ j and any geodesic arc y{t), < t < 1, 
connecting Xi to Xj we obtain: 

|(/)(xi) - (l){xj)\ < / \^x{t)(t>\ ■ \x(t)\dt < V2DE{xi,Xj) . 
Jo 

In addition, for any £ M. 

k rk 



H{(t>;ii)<- \V<pfdn+ Vdn<E 
^ Jr 
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We now show the existence of (j) satisfying (|4.2jl . provided 

\(t>i - 4'j 



max „ , — < V2 (4.3) 



ij^j DE{xi,Xj) 

Given e > 0, let D^^{x,Xi) := max{L'£;(x, Xj), e}. Set 

(l)%x) = min I \/2i:)|;(j;, j;,) + 0j - el . 

l<i<n L J 

Note that D^^{x^Xi) is a Lipschitz function on (recall > and satisfies \\/ DE{-,Xi)\ < 
\/E — V for almost any x G M'^. In addition, the condition H4.3() implies that (j)'^{x) = if 
|x — < e for any j E {1, . . . n}, if e is sufficiently small. Now, let us take the smoothing 
kernel 

rj G C°°(M'=) ; > 0, ; /" 77 = 1 , ; 7/(2;) = |x| > 1 . (4.4) 
Let ?75(x) = 5~^rj{x/5). Let 

(^(x) := ris*4f . 

Evidently, the inequality \S/(t)'^\ < \pl\jE — V holds now everywhere for (p. If (5 < e then, by 
the last condition in 1)4. 4(1 . also the condition (j)'^{xi) = cpi are preserved for (p. So 

hC$)< inL|g,max < V^X (4.5) 

E>V I »7^i DE{Xi,Xj) J 

is verified. 

To prove the opposite inequality we construct a probability measure as follows. Let E 
defined by 

V2D-^{Xm, Xk) = \(t)k - 4>m\ ; V2D^{Xi, Xj) > \<pi - (4.6) 

Assume E > V. Let /x = fj,k,m as defined in (|3.5j) . We now calculate the minimizer cp £ B-^ 
of Jjgfc |V0pd/ifcm- Note that for any G C-^(M*') we may define rj £ [0, S^^m] via r/(s) = 
(j) iQk,m{s)) , where (7^,™ is the parameterization of the geodesic arc connecting Xm to Xk (see 
paragraph preceding (|.S.3j) ). Then 

\V(l)\'^diik,m> / |??(s)lV,m(s)ds (4.7) 



where Pk,m as defined in H3.4() . Now, the minimizer on the RHS of (|4.7() subject to the 
condition r/(0) = (pm, 1] {Sk,m) = 4'k is attained for 



rjis) =4>m+ ^ 



dt 



By it follows that 



|V0|-d/...^>^ 'g^ . (4.8) 
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must holds for any (j) G B-^. 
Next, by (ESj) and (1231): 



+V = E. (4.9) 



Multiply H4.9() by pfc^m! integrate from to S'fc^m and use 1)4.8^ to obtain 



/ 



for any i;^ G S-^. This implies the reverse inequality of (|4.5jl . provided E > V. 
Finally, we observe that the choice // = (5a;Q where V{xo) = V guarantees: 

H{^)>H(^;fi) = V. 



Let now H be the Legendre transform of H: 



□ 



H*{X) := sup < -H{lp) + ^ Ai(/>i I e M U {oo} . 
"^eK" I 1 J 

Similarly 

H*{X;H):= sup | //) + ^ Ai(/>i I G M U {oo} . 

Note that 6oi/i H and /x) are convex by definition. From Lemma l4. II we also obtain 

Corollary 4.1. H*{X) G M /or eac/i T G E'' provided Aj = 0. 

Proof. We first note that DE{xi,Xj) w \/^ as -E — > oo. From Lemma IXD it follows that H is 
at least quadratic with respect to \(pi — (pj\ for any i / j. So, if we fix, say, (pi = 0, then 

lim — =f — = oo . 

I (p |^oo,</)i=0 I <p I 

In addition, H is invariant under the shift (f> — > (f> + t 1 where 1 = (1, . . . 1) G and t G M. 
This implies that Aj^j — H{ (()) is bounded from above on M'^, provided Aj = 0. □ 

From Proposition 123 and (|2.2j) we observe that 

L{ii)=H*(X-ii) . (4.10) 
By definition, H*{ A ; /i) > H* { A ) for any ji G M.. So, we obtain: 
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Corollary 4.2. The minimal action is not smaller than H { \ ). 
We now state 

Lemma 4.2. Let (p G d^H*{-; fx). Assume H{(();fi) = H{(j)) and H{-;fj,) is a convex 
function in the first variable. Then /x is an action minimizer (that is, a minimizer of \4-i(^ )- 

Proof. Consider the chain of inequahties: 

H*(X; fi) + H{^) > H*CX; fi) + /u) > ^ • X , 

which holds for any fi G A4, (f> G M" and A in the domain of H* . The left inequality follows 
from the definition of i/ as a maximizer of H over A4, while the right inequality follows by 
the convexity of H{-; fi) and Lemma 13.31 -iii. The conditions of the Lemma and Lemma l3.3l -iii 
again imply that both inequalities are, in fact, equalities , so 



Now, by definition 



so 



H*{\-ii) + H{(t))= 4) ■ A 



Now, H < H hy definition, so 



i/*(A) + if**(0) < cP ■ X 



(4.11) 
(4.12) 

(4.13) 



holds as well. Now, Lemma l3.3t -iii implies that the reverse inequality must hold in 1)4. 13(1 . 
Hence there must be an equality in (|4.13)) . which induces the equalities in (|4.12|1 and (|4.1H1 
as well. It follows that /i is an action minimizer by (|4.1U|) and Corollarv 14.21 



□ 



Corollary 4.3. H*{X) = sup^>y {\/2wiJ^ ( T) - . 
Proof. By definition of H* and Lemma l4.H 



H (X 



mm 



min in^ 

~6 E>V 



H{<p)- X ■ 

'V2W^^\'X)-E 
where we used the duality relation given by 1)3.61) . 



E — X ■ <p ; max r 



< V2 



sup 

E>y 



(4.14) 
□ 



Lemma 4.3. Suppose Eq > V is the minimizer of \4.14^ . Then there exists A^^^ G Q-^ 

which minimize the Wasserstein cost W^^\ i.e. ^f°j^Eo{xi.,Xj) = W^^^(A), 

and 
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IS 



satisfied. If Eq = V then for any such G Q-» , the inequality 



holds. 

Proof. Let En \ Eq. For each n, set A^" S be a minimizer of the Monge-Kantorovich 

problem associated with W^^(A). Note that such a minimizer may not be unique. We 
choose a subsequence so that the hmit 

A^o ■- lim A^" (4.15) 

n— »oo 

exists. Evidently, A^o is a minimizer of the Monge-Kantorovich problem associated with 
W^^|( A ) (again, possibly not the only one). 
Next, since Eq is a maximizer of ()4.14|) . 



V2 (W«( A ) - Wg( A )) <En-Eo 



so 



E E ^& iDEAx^,X,) - De,{x,,X,)) < E E ^f]DEM,X,)-Y, E ^^o(^.,^ 

iei+jei- i€i+jei- i&i+j&i- 



= W^^^( A ) - W^;^( A ) < . (4.16) 

Take the limit n — > oo and use 1)4. 15() and Lemma l3. II to obtain from 1)4. 16() 



E Y.KjT^AEo)<V2. (4.17) 

iei+jei- 

This completes the proof for the case Eq = V. 

Now, let En Z' Eq, and let us consider the subsequence along which the limit 

A^o ■- lim A^" (4.18) 

Jl— >oo 

exists. By following the preceding argument we obtain 

E E <^ ^ ^ • (4.19) 

Finally, if the maximizer A^'^ of W^-*^^ ( A ) is unique, then Af° = Af° := A^o and the 

proof follows from 1)4.171 UnjU]) . Otherwise, since both A-^ are minimizers, so is the convex 
combination thereof. Now, we utilize (|4.17l 14.19(1 to choose a G [0, 1] for which A^o := 
0(A-^ + {1 — ce)A-^ satisfies the desired equality 

wg(^) -EE <F^AEo) = V2 . 
i&i+jei- 

□ 
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Given G M" so that H{4>) = E. Let 

Ie ■■= {(i,i) el+xl^ ■ - = V2DE{xi,Xj)^ . (4.20) 

Let Ce be the convex hall of the set of measures /ijj as defined in ()3.5() where G Ie, 
that is 

fi£ Ce iff a* = X] ; "ij > ' X] X] ^ ■ i'^-'^'^) 

iei+jei^ iei+j&i- 

Lemma 4.4. For any </> G satisfying H{ (f>) = E > V 

h(1) = h{^;Ii)=E 

holds for any fj. £ Ce- 

Proof. The proof of Lemma 14.41 is, basically, identical to the proof of the second part of 
Lemma l4. II □ 

Lemma 4.5. Let E > V and /x G Ce for some cj) G M" satisfying H{4>) = H{(f>;^) = E. 
Then for any C 

In particular, II{-,^) is convex for any fx of the form \4.21\ ). 

Proof. Given C G M", let C £ C^(M'^) be an optimal solution corresponding to /i. We push 
it backward to a function on the graph composed of Uij[0, Sij] via 

Vidis) := ({qij{s)) , 

so ^ 

H{C; /") < J 5" "ij / Pi,i(5)l'?ijf + [vdfi. (4.22) 

The equality in 1)4.22(1 is achieved if we minimize the RHS on each branch separately, subjected 
to the prescribed end conditions r]ij{0) = VijiSij) = Cj- Hence Pijrjij is constant on each 
branch. Taking the end conditions and the definition of pij ()3.4I) we obtain that 

hM) - {£■' ^) " (0 - C.) - 5^^!^ . (^.23) 

Now, we perturb C« ~^ Ci + ^- Let be defined on the {i,j) branch [0, Sij] where C(0) = 1, 
C{Sij) = 0. Let Ci be the canonical i— unit vector in R", then, evaluating H( ( +eej; /u) via 
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the function r/j.j + eQj on each branch [0, Si,j] where aij > (fixed i) and integration by 
parts yields 

HiC +eef,fj.)-H{C;fi)=eY,ai,jPi,j{0)iiiji0) + O{e'^) , 

j 

From (|^^ it follows that 

and the proof follows by integration. □ 
Corollary 4.4. The function H is convex on M.^. 

Proof. By definition, H{ (j) ) is the maximum of the set H{ cj) ; /i) where run on M.. From 
Lemma 14.51 it follows that this maximum is obtained at a convex function. Hence, H is the 
maximum of a family of convex functions, so it is convex. □ 



Proof of the Main Result: 



We prove the main result in full generality (y ^ 0). 
We obtain from (|4.10j) and Corollarv 14.31 that (|1.21j) is a lower bound for the minimal action. 

By Corollarv 14. 11 A is in the domain of H* if Aj = 0, so there exists (j) £ ^"a 
H{ (f> ) = E and assume that E > V . By Lemma l4.ll 



max 



DE{xi,Xj] 



< V2 . (4.24) 



Since H*{ X) + E = H{ 4>) + H*{ A ) = ^ Aj^^ we obtain the equality in (|4.25|) below from 
Corollarv Ol 

n 

Wg)(t) = 2-1/2^^^^. . (4.25) 

1 

Let 

/^ = 2-'^'E E A':^.^'.:'(^Ki (4.26) 
iei+jei- 

where fiij as defined in (|3.41 13. 5|) . corresponding to the same energy E. By Lemma ESI is 
a convex combination of fXij, hence a probability measure. 
Now let A G Q^. By (irailOKll 

^ A,,jDE{xi,Xj) > 2-1/2 ^ ^ ^. .|^^ _ ^^.| 
iei+jei- iei+jei- 

n 

i£l+jel^ 1 
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Now substitute the minimizer G Q-r» for A above. Then the chain of inequahties turn 
into equahties. In particular we obtain that, for any i € I+, j G I~, either Afj = or 
(pi — = \/2DE{xi, Xj). As a result we can apply Lemma 14.41 to obtain 

H{^;fi) = E. (4.27) 

In addition, we substitute the equalities — (f)j = \^DE{xi, xj) and aij = 2^^^'^AfjTij{E) 
in Lemma 14.51 to obtain 

V^i7(-,/.) = T, 

which implies 

7 ed^H*{-,fi) (4.28) 

By (|4.27l4.28j) and Lemma 14.21 we obtain that is an action minimizer. In particular, it 
follows that (|1.2H) is the minimal action. 

Finally, if E = V, then let (3 = 2-^/^Y.iei+ ^fj'^idi^)- % lemma Ol /3 < 1. 
define 

= 2-'/' E E + (1 - P)K (4.29) 

ie/+ j€l- 

where xq is a maximizer of y, so V{xo) = V. Since E = V the equality (|4.27|) holds for /u 
given by 1)4. 29() . In addition, 1)4. 28() is also verified for this /x by Lemma 14.51 Hence, /x is an 
action minimizer via Lemma 14.21 as well. 



5 Conclusion 

We considered the extended minimal action principle for stationary actions in the presence of 
point sources and sinks. The main conclusion of this paper is that this minimal, stationary 
action is obtained as a minimization of a metric Monge-Kantorovich for the (non-normalized) 
pair of discrete measures ()1.14() . This stands in contrast to the non-metric Monge-Kantorovich 
transport (|1.5|1.6|1 which is related to the extended minimal action in the non-stationary case. 

Another interesting conclusion is the relation H1.23() between the expectation of the 
inverse-time of the minimal orbits to the given flux. This equality follows from substitu- 
tion of (|1.22|) and using the fact A^ G in the middle term of (|1.23|) . We stress that this 
relation is preserved in the case E = V, since the point measure PS(^xo) corresponds to orbits 
of infinite time length-these are the orbits which converges to the maximum of V, but never 
get there, since the metric \/E — V ds is degenerate at the point xq where V{xq) = V. 
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